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A consistent set of thermochemical property values, A/#°, A/G°, S", and C' p , at 298.15 K is given for the 
known constituents of aqueous sulfur dioxide (SO§(aq), HSOj(aq), SO^-<aq), H + (aq), and S 2 0?"(aq)). Also 
tabulated are values of the mean ionic activity coefficients, osmotic coefficients, partial pressure of S0 2 (g), and 
the relative apparent molar enthalpy as a function of concentration of SOj(aq) at 298.15 K. The data analysis 
considered a wide variety of measurement techniques: calorimetric enthalpies of solution and reaction, heat 
capacities, equilibrium constants, solubilities, and vapor pressure measurements, both partial and total, over 
aqueous solutions of SOi for the temperature range 278 to 393 K. All auxiliary data have been taken from the 
most recent set of COD ATA values which were converted to a standard state pressure of one bar (0.1 MPa). 
For the process S0 2 (g) = SO§(aq), the selected "best" values are: if = 1.23 ±0.05 mol kg" 1 bar" 1 , 
AG°= -0.51 ±0.10 kJ mol-', Aif=~26.97±0.30 kJmol- 1 , and AC;=155± 10 J mol" 1 K- 1 . The standard 
state partial molar entropy of SOj~(aq), obtained by the analysis of data via two independent thermodynamic 
pathways is -15.40±0.8O J mol- 1 K"' at 298.15 K. Parameters are given which extend the predictions to 
temperatures up to 373 K. 
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1. Introduction 

The purpose of this paper is to present the results of an 
analysis of the thermodynamic properties of the 
S0 2 + H 2 system. This study was performed as part of 
the analysis of the thermodynamics of sulfur and its 
compounds. The aqueous species of relevance to this 
study are neutral aqueous sulfur dioxide SO^aq), bisul- 
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fite ion HSOj(aq), sulfite ion SOi"(aq), pyrosulfite ion 
S 2 Os~(aq), and H + (aq). A consistent set of property 
values 1 , A f H°, A f G°, S°, and C° at 298.15, is given for 
these species which results from the evaluation of the 
various processes involving them. Also included are 
values of the activity and osmotic coefficients, the rela- 
tive apparent molar enthalpy (£4,) as a function of con- 
centration, and the partial pressure of S0 2 (g) over its 
aqueous solution. 

Throughout this paper we are adhering to the stan- 
dard states and units used in the NBS Tables of Chem- 
ical Thermodynamic Properties [l] 2 . All auxiliary ther- 
mochemical data have been taken from the most recent 
set of CODATA key values [2] which were converted, 



' See Glossary of symbols at end of paper. 
2 Figures in brackets indicate literature references. 
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where appropriate, to a standard state pressure of one 
bar (0,1 MPa). 

Thermodynamic data from a wide variety of meas- 
urement techniques are considered. They include cal- 
orimetric measurements of heat capacities and en- 
thalpies of solution and reaction, measured equilibrium 
constants, entropies, solubilities, and vapor pressure 
measurements (both partial and total) over aqueous 
solutions. The approach adopted in this paper is to select 
a tentative set of values and accompanying uncertainties 
for the various processes involving sulfur (IV)-oxygen 
species and then to examine the various thermodynamic 
pathways for consistency. A final set of values for both 
the processes and thermodynamic properties is then se- 
lected which provides a "best" fit to all of the available 
data. The final uncertainties are adjusted to reflect the 
agreement or lack of it between the various data sets. 



Na 2 S0 3 (cr) = 2 Na + (aq)H-SOi"(aq) 



(G) 



together with the Third Law entropy for Na 2 S0 3 (cr), 
the enthalpy of solution, and the value of AG for the 
solution process (G). 



3, Equilibrium Calculations on Aqueous 
Solutions of Sulfur Dioxide 

The equilibrium constants for the equilibria in aque- 
ous solutions of S0 2 as described in processes (A), (B), 
and (C) are: 



K A =[a{H + (zq)} 5{HS07(aq)}] 
^[a{SO?(aq)}a{H 2 0(/)}] 



(1) 



2. Key Processes Involved in the Evaluation 

The species considered to be present in an aqueous 
solution of S0 2 (stoichiometric "S0 2 (aq)") are neutral, 
unionized SO§(aq), HSOl(aq), SO! _ (aq), S 2 0?"(aq), and 
H + (aq). They are involved in the following equilibria: 

SO?(aq) 4- H 2 O(0 = H + (aq) + HSO 7 (aq) (A) 



tf B =S{H + (aq)} a{SO?-(aq)}/5{HS07(aq)} (2) 



HSOy(aq) = H + (aq) + SO§-(aq) 



2 HSOT(aq) = S 2 Ol-(aq) + H 2 0(/) 



(B) 



(Q 



Selections are made for the various parameters that 
enter into the solution process of S0 2 (g) into water so 
that a set of "best" values of AG and AH can be 
obtained for the processes: 

S0 2 (g) = SO?(aq) (D) 

S0 2 (gHH 2 0(/)=H + (aq)+HS07(aq) (E) 

The standard state partial molar entropy of SO^(aq) 
can be obtained from the standard state Gibbs energy 
and enthalpy changes for the process: 



S0 2 (g) + H 2 0(/) = 2 H + (aq) + SOl-(aq) 



(F) 



which is the summation of processes (A), (B), and (D). 
This partial molar entropy can also be obtained indepen- 
dently using the data for the process: 



and 



tfc=a{S 2 Oi-(aq)} 5{H 2 0(/)} 
-h3 2 {HSOT(aq)}. 



(3) 



In the above equations, 3, is the activity of the i" 1 species 
which is equal to the quantity rhj^j where fh\ and y are, 
respectively, the molality and activity coefficient of that 
species. The symbol " " " denotes a quantity which 
pertains to a species as distinct from a stoichiometric 
quantity [3], The absence of the symbol " " " over a 
quantity implies that the quantity is stoichiometric. The 
treatment of the experimental data often requires some 
assumptions about the activity coefficients of the species 
in solution and the mathematical solution of these three 
simultaneous, nonlinear equations for the molalities of 
all of the species considered to be present in solution. 
The first requirement will now be addressed. 

Since we have no direct knowledge of the activity 
coefficients of these or any other individual aqueous 
charged species, it is necessary to make assumptions to 
obtain values for these activity coefficients. It is assumed 
that the activity coefficients can be calculated using the 
expression: 



lny i ^-A m I m ^/(l+BI wl ) 



(4) 



where A„ is the Debye-Hiickel constant, B is an "ion- 
size" parameter, and I is the ionic strength which is 
calculated as: 



/=(l/2)2/M. 



(5) 
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The above summation extends over all of the species in 
the solution. The Debye-Hiickel constants tabulated by 
Clarke and Glew [4] were used in our calculations. Sta- 
ples and Beyer [5] have recently calculated values of the 
Debye-Hiickel limiting slopes using the IAPS recom- 
mended values for the dielectric constants and density 
of water as a function of temperature [6,7]. These new 
limiting slopes differ by 0.24 and 0.61%, respectively, 
from the Gibbs energy and enthalpy limiting slopes 
given by Clarke and Glew [4]. The effect on our final 
recommended values due to uncertainties in the Debye- 
Hiickel slopes is negligible in comparison with the ex- 
perimental errors in measured properties which were 
extrapolated to zero ionic strength whenever possible. 

The B parameter in eq [4] is a quantity which serves 
to account for the repulsive forces between the par- 
ticles. If sufficiently accurate activity or osmotic coeffi- 
cients were available for aqueous solutions of S0 2 , it 
would be possible to regress a value for this parameter 
using a chemical equilibrium model [3]. Since this infor- 
mation is not available, it is necessary to work with 
estimated or inferred values of B. Evidence for a value 
of B equal to 1.5 was obtained when extrapolating mea- 
sured values of K A (sec. 4.1) to zero ionic strength; spe- 
cifically, the extrapolation could frequently be made 
with a line of lesser slope when a value of 1,5 for B was 
used to calculate values of %. Also, a B value of 1.5 is 
typical of the values obtained by Hamer and Wu [8] 
from their correlation of the activity and osmotic coeffi- 
cient data of aqueous uni-univalent electrolytes. Evi- 
dence for a value of B equal to 2.5 was found in the 
treatment of the heat of solution and heat of dilution 
data (sec. 4.2) to obtain values of LH&. There the best fit 
of the experimental data was frequently obtained using 
a value of B equal to 2,5. Less direct evidence for a value 
of 2.5 comes from a treatment of the osmotic coefficients 
of aqueous sulfuric acid using an equilibrium model [3] 
in which it was found that a value of B equal to 2.5 could 
represent the osmotic coefficients of sulfuric acid to 
within 1.2% up to a stoichiometric molality of 0.20 mol 
kg~~'. This stoichiometric molality corresponds to a cal- 
culated (see eq (5)) (species) ionic strength of 0.25 mol 
kg -1 . Since an aqueous solution of S0 2 with a stoichio- 
metric molality of 2.0 mol kg -1 has a calculated ionic 
strength of only 0.21 mol kg"', the use of assumed val- 
ues of B ranging from 1,5 to 2.5 will allow us to proceed 
with the necessary equilibrium calculations without 
serious error. 

The solution of the simultaneous, nonlinear equations 
(1), (2), and (3) was accomplished using an iterative 
numerical procedure (C05NCF) which is a part of a 
library of numerical analysis programs [9]. The calcu- 
lation was made self-consistent both in regards to the 
activity of the water which is a participant in the equi- 



ilibria (see reference [3]) and with eqs (4,5). This model 
allows us to calculate the amounts of each species in the 
solutions and values of the fractions (ai) of the species, 
i.e., a, = hi/m st where m st is the stoichiometric amount of 
S0 2 in solution calculated as m s '(H2S0 3 ) for con- 
venience. It does not imply the physical existence of 
H 2 S0 3 . Results of these calculations at 298.15 K are 
given in table 1 using our final selected values of K K , K B , 
and K c and a value of B equal to 2.0. 

Throughout, where necessary and possible, reported 
equilibrium constants have been corrected to 298.15 K, 
for the activity, and to zero ionic strength. 

4. The Thermodynamic Parameters for 

the Description of the Equilibria 

in Solution at 298.15 K and as a 

Function of Temperature 

4.1 Process (A) 

The reported equilibrium constants at 298.15 K for 
process (A), the first ionization of SO°(aq), are sum- 
marized in table 2a. The most detailed investigations 
cited are those of Tartar and Garretson [15] and Huss 
and Eckert [24]. While previous critical evaluations 
[1,32] have relied upon the results of Tartar and Gar- 
retson [15], we have selected the data of Huss and 
Eckert [24], ^T A =0.0139±0.004 a mol kg -1 , as being pref- 
erable since it is based upon two different measurement 
techniques (conductivity and ultraviolet spec- 
trophotometric measurements) which avoid the possible 
systematic errors inherent in the use of electrochemical 
cells which were used by Tartar and Garretson [15]. 

Measurements of K A as a function of temperature 
have been used to obtain values of AH a which are cor- 
rected to 298.15 K (see table 1) using the model of 
Clarke and Glew [33]. A value of AC; = -272±10 J 
mol -1 K _1 for process (A) was used in performing these 
calculations. This heat capacity change is based upon 
the calorimetrically determined standard state heat ca- 
pacities of SO?(aq) and HSOj(aq) of, respectively, 
195± 10 J mol" 1 K-' (Barbero et al. [34]) and -2± 10 J 
mol" 1 K _1 (Allred et al. [35]). The uncertainties as- 
signed to the calculated enthalpies in table lb are purely 
statistical and refer to 95% confidence limits. 

The model of Clarke and Glew [33] is based upon a 
Taylor series expansion of the heat capacity at a refer- 
ence temperature (298.15 K). When dAC°/dT is con- 
stant over the temperature range of interest, the 



* The thermodynamic equilibrium constant is dimensionless. The units expres- 
sion identifies the composition units used for components in K. 
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Table 1. Fractions (a) of species and the species ionic strength (?) at 298.15 K in aqueous sulfur dioxide solutions as a function of the 
stoichiometric molality of S0 2 (>m s, (S0 2 )) and of (hypothetical) H 2 S03(m 5l (HS0 5 )). These values were calculated using the final selected 
values of K At K B , and K c and a value of B equal to 2.0. 



m s, (S0 2 ) 


m"(H 2 S0 3 ) 


a(H+) 


a(HSOj) 


a(SOi-) 


a(SOS 


2a(S 1 Ol-) 


/ 


mol kg 


-i 












mol kg-' 


0.001000 


0.001000 


0.9406 


0.9404 


0.7697 XlO" 4 


0.0595 


0.6058 XlO- 4 


0.9407X10- 3 


0.002000 


0.002000 


0.8949 


0.8947 


0.4036X10- 4 


0.1051 


0.U23X10" 3 


0.1790X10- 2 


0.003000 


0.003000 


0.8576 


0.8574 


0.2781 XlO" 4 


0.1424 


0.1 573 XlO" 3 


0.2573X10- J 


0.004000 


0.004000 


0.8261 


0.8259 


0.2141 XlO" 4 


0.1739 


0.1971XL0- 1 


0.3305 XlO" 2 


O.OOSOOO 


0.005000 


0.7989 


0.7986 


0.1749X10" 4 


0.2011 


0.2329X1O" 3 


0.3995 XlO" 2 


0.006000 


0.006001 


0.7750 


0.7747 


0.1485X10- 4 


0.2250 


0.2654X 10- 3 


0.4651 XlO- 2 


0.007000 


0.00700 i 


0.7537 


0.7534 


0.1293X10" 4 


0.2463 


0.2952X10" 3 


0.5277 XlO- 2 


O.OOSOOO 


0.008001 


0.7345 


0.7342 


0.1 148X10" 4 


0.2655 


0.3227X1O- 3 


0.5878 XlO" 2 


0.009000 


0.009001 


0.7172 


0.7168 


0.1034X10" 4 


0.2828 


0.3483X10- 3 


0.6457X10- 2 


0.01000 


0.01000 


0.7013 


0.7009 


0.9412X10" S 


0,2987 


0.3723 XlO' 3 


0.7016X10" 2 


0.02000 


0.02000 


0.5919 


0.5913 


0.5111X10-' 


0.4081 


0.5524X10" 3 


0.1185X1O-' 


0.03000 


0.03001 


0.5271 


0.5264 


0.3589X10"' 


0.4729 


0.6742X1O" 3 


0.1583X10-' 


0.04000 


0.04002 


0.4823 


0.4816 


0.2797X10- 5 


0.5177 


0.767OX 10" 3 


0.1932X1O-' 


0.050O0 


0.05004 


0.4487 


0.4479 


0.2307X10- 5 


0.5513 


0.8422X10- 3 


0.2247X10"' 


0.06000 


0.06006 


0.4221 


0.4212 


0.1972X10" 5 


0.5779 


0.9056X10- 3 


0.2537X10"' 


0.07000 


0.O7O08 


0.4003 


0.3994 


0.1728x10"' 


0.5997 


0.9605 XlO- 3 


0.2807X10"' 


0.08000 


0.08011 


0.3821 


0.3810 


0.1541 XlO" 5 


0.6179 


0.1009 XlO" 2 


0.3062 x 10"' 


0.0900O 


0.09014 


0.3664 


0.3653 


0.1393 X10" 5 


0.6336 


0.1052X1O" 1 


0.3304x10"' 


0.1000 


0.1001 


0.3527 


0.3516 


0.1273 XlO" 5 


0.6473 


0.1092 XlO" 2 


0.3535X10"' 


0.1500 


0.1504 


0.3035 


0.3022 


0.9015X10-" 


0.6965 


0.1248X10- 2 


0.4565X10" 1 


0.2000 


0.2007 


0.2718 


0.2705 


0.7065 XlO" 6 


0.7282 


0.1363X10" 3 


0.5456X10"' 


0.3000 


0.3016 


0.2318 


0.2303 


0.5O15X10" 6 


0.7682 


0.1532X10- 2 


0.6987 x 10"' 


0.4000 


0.4029 


0.2066 


0.2049 


0.3934X10" 6 


0.7934 


0.1655X10- J 


0.8308X10"' 


0.5000 


0.5045 


0.1887 


0.1869 


0.3260X10- 6 


0.8113 


0.1754X1O" 2 


0.9493x10"' 


0.6000 


0.6065 


0.175! 


0.1732 


0.2796X10" 6 


0.8249 


0.1836X10- 2 


0.1058 


0.7000 


0.7089 


0.1643 


0.1624 


0.2457 XlO" 6 


0.8357 


0.1 906 XlO- 2 


0.1159 


0.8000 


0.8116 


0.1554 


0.1534 


0.2195X10" 6 


0.8446 


0.1968X10" 2 


0.1254 


0.9000 


0.9148 


0.1480 


0.1459 


0.1988X10-* 


0.8520 


0.2023 XlO' 2 


0.1344 


1.0000 


1.0183 


0.1416 


0.1395 


0.1819X10" 6 


0.8584 


0.2072 XlO' 2 


0.1430 



temperature dependency of the equilibrium constant is 
given by: 

Rln K= -AG 2 98 .ls/298.15 + Aff 3 V» s (l/298.15-l/7 , ) 
+ AC; [(298. 1 S/D - 1 + In (J/298. 1 5)] 
+(298. 15/2) (dAC;/dr)[(77298.15)- 
(298. 1 5/jT)- 2 /n (77298,15)] (6) 

A convenient aspect of the above equation is that known 
values of AG , MI", AC°, and (dAC p °/dr) at a reference 
temperature provide all of the information needed to 
calculate the variation of the equilibrium constant with 
temperature. 

Given in table 2b are enthalpy values obtained from 
the direct calorimetric measurements of Vanderzee [41], 
Dobrogowska and Hepler [39], and Zambonin and Jor- 
dan [38] and from the enthalpies of solution of S0 2 (g) as 
a function of concentration from the measurements of 
Johnson and Sunner [37] and Stiles and Felsing [36]. 
Values of A/f£ were also obtained from the heat of 
dilution measurements of Dobrogowska and Hepler [40] 



using a regression calculation in which AHl was varied 
in the equilibrium modeling calculation of L% (sec. 6) 
until agreement was obtained between measured and 
calculated values of the heat of dilution. A similar pro- 
cedure was used for treating the aforementioned en- 
thalpy of solution measurements of S0 2 (aq) where val- 
ues of L$ could be obtained from the concentration 
dependence of the heat of solution measurements. We 
adopt a value of AHI= - 17.80±0.40 kJ mol~' based 
largely on the direct calorimetric measurements of Dob- 
rogowska and Hepler [39] and ofVanderzee [41] and the 
enthalpy of solution measurements of Johnson and 
Sunner [37]. Note that in table lb the direct calorimetric 
measurements of Dobrogowska and Hepler [39] and of 
Vanderzee [41] are sensitive to the B parameter. This 
causes an uncertainty in the final value of MIX compara- 
ble to the errors in the measurements themselves. 

4.2 Process (B) 

The thermodynamic data for process (B), the second 
ionization for SO^aq), are summarized in table 3. 
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Table 2. Thermodynamic parameters for process (A): SO?(aq)4- H 2 0(/) = H + (aq)4-HSO7(aq). 



a. Equilibrium constant at 298.15 



b. Enthalpy change at 298.1S K 



Investigator(s) 


tf/mol kg- 1 


Investigator(s) 


A#°AJ mor 1 


Drucker [10] 


~0.009 ,b 


Stiles and Falsing [36] 


-17.2±0.7' 


Sherill and Noyes [11] 


0.01 2 C 


Johnstone and Leppla [13] 


-16.4+0.9 


Britton and Robinson [12] 


~0.008 a 


Johnson and Sunner [37] 


-18.0±0.6 a 


Johnstone and Leppla [13] 


0.01 30 d 


Deveze and Rumpf [19] 


-16.7±1.0 


Yui [H] 


0.0127 


Flis et al. [20] 


-16.7±10. 


Tartar and Garretson [15] 


0.0172 


Zambonin and Jordan [38] 


-14.3±1.5 


Frydman et al. [16] 


0.02 3 


Dobrogowska and Hepler [39] 


-17.4±0.3 b 


Ellis and Anderson [17] 


0.014 


Dobrogowska and Hepler [40] 


-L8.2±1.0 a 


Ludemann and Franck [18] 


0.0103 e 


Vanderzee [+1] 


-17.8 2 ±0.3 C 


Deveze and Rumpf [19] 


0.0145 






Flis et al. [20] 


0.010 






Deveze [21] 


0.013 






Sekine et al. [22] 


0.014* 






Beilke and Lamb [23] 


0.0165" 






Huss and Eckert [24] 


0.0139 







Method 



calorimetric-heats of solution 

KvsT 

calorimetric-heats of solution 

KvsT 

KvsT 

calorimetric-heats of solution 

calorimetric-neutratization 

calorimetric-heats of dilution 

calorimetric-neutralization 



* Corrections were applied. 

b Drucker's [10] results are based on the mea- 
surements of McRae and Wilson [25] and of 
Walden and Centnerzwer [26]. 

* Based on measurements of Kerp and Bauer 
[27] and of Lindner [28]. 

a Based on measurements of Campbell and 
Maass [29] which are also reported in Beazley et 
al. [30] and Morgan and Maass [31]. 

" Pressure is 0.2 kbar (20 MPa). 



"These values were obtained by regression calculations. Additional information obtained from 
these regressions are as follows: Stiles and Felsing [36], (21 points) B=2.5, LH$= -27.14±0.10 kJ 
mol- 1 . Johnson and Sunner [37], (7 points) B = 2.5, A^=-26.W±0,04 kJ mol" 1 . Dobrogowska 
and Hepler [40], (21 points) B=2.0 to 2.5. The uncertainties given here refer to two standard 
deviations as distinct from the overall assigned uncertainties given above. 

b This value was obtained using the Guggenheim equation with ji= —0.2. If an extended Debye- 
Huckel equation is used with B = 2. 5, a value of —17.17 kJmol -1 is obtained; if B is set equal to t.5 
then 4//°=- 17.61 kJmo]- 1 . 

c This value was obtained using an extended Debye-Hiickel equation with B= 1.6. If a value of 
B = 3.0 is used, a value of &H£ equal to — 17.43 kJ mol" ' is obtained. 



Table 3, Thermodynamic parameters for process (B): HSOy(aq) = H + (aq) + S03~(aq). 



a. Equilibrium constant at 298,15 K 


b. Enthalpy change at 298.15 K 




Investigator(s) 


^XlOVmolkg" 1 


Investigator(s) 


AtfVkJ mol" 1 


Method 


Jellinek [42] 


500' 


Arkhipova et al. [45] 


-12.1 ±4.0 


KvsT 


Kclthoff [43] 


~2 b 


Zambonin and Jordan [S&] 


-5.2±0.8 


calorimetric 


Britton and Robinson [12] 


~4" 


Hayon et al. [46] 


-11.6+3.3 


KvsT 


Yui [14] 


6.3 


Teder [47] 


-8.0±5.0 


KvsT 


Tartar and Garretson [15] 


6.24 


Krunchak et al. [48] 


-16.7±6.7 


KvsT 


Frydman et al. [16] 


~10 b 


Allred et al. [35] 


-3.6±0.2 


calorimetric 


Cuta et al. [44] 


6.4 b 


Vanderzee [41] 


^3.67±0.07 


calorimetric 


Arkhipova et al. [45] 


5.05 








Hayon et al. [46] 


6.3 








Teder [47] 


7.9 








Krunchak et al. [48] 


6.61 









* Uncorrected. 

b Corrections were applied. The Cuta et al. [44] original extrapolation to zero ionic strength !ed to K=(7.1±0.5)X t0~ 8 mol kg~ 



From table 3a we select a value of K B equal to 
(6.5±0.5)x 10 s mol kg~'. As was done for process (A), 
values of AHb were calculated from the variation of K B 
with temperature using a constant value of 

AC° = ~262±14 



SO!-(aq) = 



J mor 1 KT 1 (c; of 

264±10 J mol- 1 K- 1 from Allred et al. 
[35]) and the Clarke and Glew model. The values of 
AHg obtained in this way are in sharp disagreement with 
the recent calorimetric measurements of Allred et al. 
[35], Vanderzee [41], and Zambonin and Jordan [38], 



The calorimetric values are preferred and a value of 
-3.65+0.10 kJ mol" 1 is adopted for LH% based on the 
first two calorimetric values. Thus for process (B) at 
298.15 K, the tentative selected values are: 
A'=(6.5±0.5)Xl0- 8 mol kg -1 . AG°=41.02±0.20 kJ 
mor 1 , AH =-3.65 ±0.10 kJ mor 1 , and AC£ = 
— 262±14 J mor 1 K _1 . It will later be necessary to 
modify AG to 40.94 ±0.20 kJ mol -1 to obtain thermo- 
dynamic consistency in a thermodynamic network, re- 
sulting in a predicted K B of (6.72±0.5)X 10 -8 mol kg -1 . 
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4.3 Process (C) 

The relatively few available equilibrium data for pro- 
cess (C), the formation of S 2 Os _ (aq) from HSOT(aq), are 
summarized in table 4a. The molar absorbance of 
S 2 05 _ (aq) obtained by Connick et al. [52] was used to 
calculate values of K c from the measurement data ob- 
tained by the earlier workers [49-51]. A value of 
K c = 0.032+0.01 moL 1 kg is adopted from the mea- 
surements of Connick et al. [52] who also obtained 
A77£=4.6±4.0 kJ mol" 1 (table 4). A value of 
AC° = -21+32 J mol-' K" 1 is estimated for process 
(C) based upon an estimate of - 100 ±25 J mol -1 K~' for 
the partial molar heat capacity of S 2 Os~(aq). Thus for 
process (C) at 298.15 K: ^T=0.032±0.0I mol - ' kg, 
AG°=8.53±0.80 kJ mol -1 , AH°= -4.6±4.0 kJ mol -1 , 
and AC P = -21±32 Jraol" 1 K" 1 . 



5. Vapor-liquid Equilibrium Data and 
the Henry's Law Constant for SO2 

There are two types of vapor-liquid-equilibrium 
measurements for the S0 2 +H 2 system which were 
considered, namely, total pressure (SO2+H2O) mea- 
surements and partial pressure (S0 2 ) measurements. 
Knowing the total stoichiometric molality of S0 2 in a 
given solution, eqs (2,3) can be solved for m{S0 2 (aq)}. 
If the partial pressure of the S0 2 (g) over that aqueous 
solution is known, the equilibrium constant for process 
(D), the Henry's Law constant, is calculated as 



K D =3{SO°(aq)}//{SC- 2 (g)} 



(7) 



In the above equation, the fugacity (f) of the S0 2 (g) 
has been taken to be equal to its pressure and the 
activity of SO°(aq) equal to its molality, i.e., 
^SO^aq^^SO^aq^c^SO^aq)}/* 8 '. The latter 
assumption is consistent with eq (4). The former assump- 
tion is justified since (i) calculations showed that the 
effect of including the fugacities as calculated from 
available second virial coefficients [53] perturbed the 



calculated values of K D by less than 1.3% up to a 
stoichiometric molality of 1.0 mol kg~\ 
(ii) the values of K D calculated from most of the data sets 
examined were extrapolated to zero mol kg -1 , and (iii) 
the scatter in the data sets does not justify this small 
correction. 

Values of the Henry's Law constant (K D ) have also 
been calculated from measurements of total vapor pres- 
sures over aqueous S0 2 solutions. To do this calculation, 
the vapor pressure of the water over these solutions was 
calculated from the activity of the water using a chem- 
ical equilibrium model [3]. The partial pressure of 
S0 2 (g) is obtained from 



P {S0 2 (g)> =p (total) -p {H 2 0(g)> 



(8) 



The Henry's Law constant was then calculated using eq 
(7). Not considered here are the results of the few in- 
vestigations summarized by references [54,55] which are 
of low precision. The Henry's Law constants obtained 
from the analysis of the data are given in table 5. Values 
of AG" and AH" for process (D) were calculated (see 
table 6) using the Clarke and Glew equation with AC° 
fixed at 155 J mol- 1 K" 1 . 

We believe that the most reliable of the vapor-liquid- 
equilibrium investigations is that of Rabe and Harris 
[62]. The reasons are: (i) they measured the partial 
pressures of S0 2 (g) rather than total pressures, (ii) the 
precision of their measurement is very good, (iii) the 
AH° calculated from their measurements is close to a 
calorimetrically determined value which will be dis- 
cussed shortly, and (iv) they took care to minimize sys- 
tematic errors due to analyses of the gas phase and aque- 
ous solutions and also due to oxidation of the solutions 
and absorption of C0 2 (g). The results of Vosolobe et al. 
[66], also based on partial pressure measurements, are 
close to those of Rabe and Harris [62]. We adopt a value 
of AG°= -0.51 ±0.10 kJ mol" 1 for process (D). This 
corresponds to a Henry's Law constant of 1.23+0.05 
mol kg - ' bar -1 at 298.15 K. The assigned uncertainty is 
large enough to overlap with the mean value obtained 
from all the AG values for process (D) summarized in 
table 6 which are based upon the partial and total pres- 
sure measurements. 



Table 4. Thermodynamic parameters for process (C): 2 HSOj{aq) = S2O^"(aq) + H20(/). 



a. Equilibrium constant at 298.15 K 




b. Enthalpy change at 298.15 K 




Investigators) K/mo\~ l kg 


Investigator(s) 


AffVkJ mol SzOHaq))- 1 


Method 


Golding [49] ~0.02 a 
Arkhipova and Chistyakova [50] -~0.02 a 
Bourneet al. [51] -~0.02 a 
Connick et al. [52] 0.032 


Hayon et al. [46] 
Connick et al. [52] 


-3.3 
-4.6 


K vs T 
K vs T 



' Corrections were applied. 
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Table5. Henry's Law constants (K) for process (D): S0 2 (g) = SO§(aq). 



T/K 


K/mo\ k<r' 


bar" 1 
[56] 


T/K K/mol kg" 1 bat- 1 


T/K K/mol kg" 1 bar- 1 
Data of Hudson [60] cont'd: 


Data of Beuschlein and Simenson 


Data of Beuschlein and Simenson [56] at 


at 0.51 g SO2/IOO g H 2 0: 






7.45 g SCV100 g H 2 0: 








308.75 




0.80 




298.35 


1.37 


321.15 


0.60 


314.15 




0.71 




304.55 


1.27 


333.15 


0.44 


320.15 




0.59 




307.35 


1.14 


343.15 


0.36 


325.15 




0.53 




310.55 


1.07 


353.15 


0.29 


335.75 




0.38 




314.15 


0.99 


363.15 


0.25 


344.15 




0.31 




317.15 


0.91 






351.75 




0.26 




320.35 


0.84 






358.15 
365.15 




0.23 
0.20 




323.55 
327.55 


0.77 
0.69 


Data of Johnstone and 


Leppla [13]: 


372.15 




0.18 




328.75 


0.66 


298.15 


1.20 


378.35 




0.16 








308.15 


0-85 


384.95 




0.14 








323.15 


0.57 


386.15 




0.13 




Data of Byerley [57]: 
















298.15 
323.15 


1.26 
0.78 


Data of Maass and Maass [61]: 


Data of Beuschlein and Simenson [56] at 
1.09 g SCV100 g H 2 0: 






283.15 


2.50 










Data of Campbell and Maass [29]: 


289.65 
295.15 


2.00 
1.63 


299.95 




1.30 




298.15 


1.20 


298.15 


1.50 


306.75 




0.98 




303.15 


1.02 


300.15 


1.38 


312.55 




0.77 




313.15 


0.75 






317.35 




0.67 




323.15 


0.57 






323.75 




0.56 




333.15 


0,45 


Data of Morgan and Maass [31]: 


334.75 




0.42 




343,15 


0.36 






340.55 




0.38 




353.15 


0.30 


273.15 


3.30 


346.55 




0.33 




363.15 


0.25 


283.15 


2.16 


352.35 




0.29 




373.15 


0.21 


291.15 


1.53 


359.55 




0.25 




383.15 


0.18 


298.15 


1.20 


371.55 
373.75 




0.22 
0.19 




393.16 


0.16 


Data oFRabe and Harris [62]: 


379.75 




0.17 




















Data of Conrad and Beuschlein [58]: 


















303.15 


1.02 










298.15 


1.17 


313.15 


0.76 


Data of Beuschlein and Simenson 


[56] at 






323.15 


0.58 


4.36 g SOj/100 g 


HjO: 










333.15 
343.15 


0.45 
0.36 










Data of Douabul and Reilly [59]: 


353.15 


0.30 


296.35 




1.28 












300.35 




1.15 




278.97 


3.40 






303.15 




1.03 




283.12 


3.03 


Data of Tokunaga [63]: 




306.15 




0.93 




288.10 


2.50 






309.15 




0.85 




292.98 


2.06 


283.15 


2.57 


311.35 




0.79 




298.15 


1.73 


293.15 


2.11 


315.35 




0.71 




303.25 


1.47 


303.15 


1.46 


318.35 




0.65 








313.15 


1.07 


321.35 




0.60 












325.15 




0.55 












328.45 




0.51 




Data of Hudson [60]: 




Data of Vosolobe et al. 


[64]: 


331.55 




0.47 












334.35 




0.44 




283.15 


2.35 


293.15 


1.40 


336.15 




0.42 




288.15 


1.85 


303.15 


1.07 


340,35 




0.38 




293.15 


1.54 


313.15 


0.80 


343.95 




0.36 




303.15 


1.07 


323.15 


0.64 


345.55 




0.35 




313.15 


0.75 


333.15 


0.53 
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Table 6. Thermodynamic parameters for process (D), SOj(g) = SO°(aq), at 298.15 K. AC£ was con- 
strained to be equal to 155 J mol~' K -1 . 



Worker(s) 


AGVkJmol- 1 


A#VkJ mol" 1 


Hudson [60] 


-0.581 ±0.029 


-27.31 ±0.29 


Maass and Maass [61] 


-O.980±0.053 


-23.5±1.9 


Campbell and Maass [29] 


-0.541 ±0.097 


-27.28±0.65 


Morgan and Maass [31] 


-0.45 ±0.1 2 


-25.6±2.2 


Conrad and Beuschlein [58] 


-0.389 




Johnstone and Leppla [13] 


-0.435±0.096 


-25.7±2.0 


Beuschlein and Simenson [56] 






at 0.51 gSO 2 /100g HjO: 


-O.66±0.15 


-30.39±0.98 


at 1.09 g SO2/IOO g H 2 0: 


-0.74±0.I2 


-28.99±0.89 


at 4.36 g SO2/IOO g H 2 0: 


-0.495 ±0.022 


-26.11 ±0.27 


at 7.45 g SOj/100 g H 2 0: 


-0.97±0.I5 


-22.6±2.4 


pooled: 


-0.72 ±0.24 


-27.0±3.4 


Rabe and Harris [62] 


-0.512±0.068 


-26.41±0.67 


Vosolobe et al. [64] 


-O.475±0.O75 


-22.2±1.2 


Tokunaga [63] 


-1.26+0.42 


-21.9±11. 


Douabul and Reilly [59] 


-1.367±0.079 


-23.8±2.0 


Byerley [57] 


-0.573 




Grand average" 


-0.57±0.10 


-25.6±1.2 



a Discarding the resulls of Douabul and Reilly [59] and of Tokunaga [63]. 



6. The Calorimetric Enthalpy of Solution 
of S0 2 (g) in Water 

Enthalpies of solution of S0 2 (g) in water have been 
measured several times [28,36,37,65-68], Roth and Ze- 
umer [68] have summarized the results of the earlier and 
not very precise investigations by Berthelot [66], Tho- 
msen [65], and Lindner [28]. These early measurements 
will not be used. Instead, only the measured enthalpies 
of solution of S0 2 (g) in water from the more recent 
investigations by Johnson and Sunner [37], Stiles and 
Felsing [36], Ramsetter and Hantke [67], and Roth and 
Zeumer [68] will be considered. The measurements 
were treated in two different ways. The first used all of 
the parameters for processes (A), (B), and (C) to calcu- 
late [3] values of the excess Gibbs energy (G ex ) as a 
function of temperature from which the excess enthalpy 
(i/ ex ) was obtained using 



H c *=G"-T(dG e VdT) p 



(9) 



The relative apparent molar enthalpy (i,^) is equal to 
H a Vm st . It is a stoichiometric quantity. 

In table lb values of AH£ were obtained from the data 
of Johnson and Sunner [37], Stiles and Felsing [36], and 
Dobrogowska and Hepler [40]. These values were ob- 
tained using a minimization, or least-squares calculation. 
As a byproduct of that calculation, values of the B pa- 
rameter were also inferred. 

Using the previously adopted value of t±Hl= — 17.80 
kJ mol" 1 and a value of B equal to 2,0 in the above 



procedure for calculating L^, the enthalpy of solu- 
tion data were treated to obtain a value of 
A#iS=-27.00±0.30 kJ mol"'. Figure 1 shows the 
graphical treatment of the data; the calculations are 
shown in table 7. To obtain thermodynamic consistency 
with later calculations, this value is adjusted to 
-26.97 ±0.30 kJ mol -1 . A value of Mfi 
= -44.77±0.50 kJ mol" 1 is obtained. 

The above procedure is essentially equivalent to using 
the experimental heat of solution data at 298.15 K and 
correcting for the enthalpies of ionization of all of the 
species. Thus 



AH aitu =AHg + [l -a{SO° z (aq)}][MfZ 
+L 4) (H + .HSOl(aq))] 



(10) 



where .L,j l (H + 'HS07) is the relative apparent molar en- 
thalpy for the (hypothetical) solution consisting only of 
the ions H + (aq) and HSOl(aq). This I.* is different than 
the stoichiometric L$ above. This type of procedure has 
been used previously by Wu and Young [69]. In eq (10) 
contributions to the enthalpy of S 2 Os~(aq) and SC>3~(aq) 
were neglected since they are negligible in comparison 
to that from HSOy(aq). In applying eq (10) below the 
electrostatic contribution [L,f,(H + -HSOl) is negligible in 
comparison to the other terms and was neglected. 

An alternative procedure is to use the above re- 
lationship (eq (10)) and to plot the experimental en- 
thalpies of solution as a function of a. In such a plot, the 
slope yields a value of AHl and the intercept at a=l 



348 



o 

to. 

CV2 

I 

m 
cd. 

CM 



o 
cm 



O 

< 



m 
f--l 

CM 

I 

O 

cd. 

CM 

I 

iO 

od. 

03 

■ l 

o 

' oi . 

CM 



LO 
OJ. 

w 

i 

o 

d. 

to 

I 

d 

CO 

I 
q 

CO 



+ Op 



a 

a 
a 



A 

a 



D 



, + fi 



# 



ff S 3 



Figure 1- Values of Aff° at 298,15 
for process (D), S0 2 (g) = 
SO?(aq), calculated from the 
heat of solution measurements of 
Johnson and Sunner [37] ( + ), 
Stiles and Felsing [36] [□], Ram- 
stetter and Hantke [67] (o), md 
Roth and Zeumer [6S](A). The 
final selected value of 
hH& is-26.97±0.30 kJ mol" 1 . 







0.0 



-! — 

0.2 



0.4 0.6 0.8 

molality/ (mol kg ~ ) 



1.0 



L2 



yields a value of Ai^g.The results of Stiles and Felsing 
[36] lead to AJ^=-17.1±1.2 kJ mol -1 and AH£ 
= — 27.1 ±0.8 kJ mol -1 . The Johnson and Sunner [37] 
measurements, neglecting the three data points at the 
lowest concentrations, lead to Ai/£= — 18.00+0.34 kJ 
mol"' and A#S= -26.76±0.40 kJ mol 1 . This indepen- 
dent computational method serves as additional con- 
firmation of the correctness of our selected values of 
A#a and AH&. The enthalpies of solution as a function 
of a(SO°(aq)) are shown in figure 2. 



7. The Properties of SO^-(aq) 

The summation of processes (B) and (E) is: 

S0 2 (g) + H 2 0(/) = 2tf + (aq) + SOr(aq) (F) 

This process is useful for obtaining the thermodynamic 
properties of SC>3 - (aq). 



7.1 The Enthalpy of Process (F) 



The enthalpy changes already obtained for processes 
(B) and (E) lead to a value of ArV£= -48.42+0.50 kJ 
mol -1 . Additional, direct experimental information 
leading to this enthalpy change is available from mea- 
surements on the enthalpy of solution of S0 2 (g) in either 
aqueous NaOH or KOH, represented as: 

S0 2 (g) + 2 0H-(aq) = SOr(aq)+H 2 0(/) (Fl) 

Use of the relative apparent molar enthalpy data of Van- 
derzee and Noll [70] for aqueous Na 2 S0 3 and the rela- 
tive apparent molar enthalpy data tabulated by Parker 
[71] for aqueous KOH and NaOH leads to values of 
-164.4+0.23, -162.3+0.30, and -161.1+2.5 kJ 
mol -1 for A#fi from the measurement data of Ram- 
stetter and Hantke [67], Roth and Zeumer [68], and 
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Table 7. Enthalpy of solution data for S0 2 (g) in water at 298.15 K. The experimental results (AH m5ils ) are given in column six and the 
corrected results, Atf ° for process (D), are given in column seven. b.H& has been constrained to equal - 17.80 kJ mol -1 . The B parameter 
was set at 2.0. 



m^HiSOj) 

mol kg' 1 



a(HSOy) 



a(SOj ) 



a(SOS) 



2 a(S 2 0?r) 



A# m 



AH" 



kJ mol 



Data of Stiles and Felsing [36] 



0.02594 

0.03839 

0.04074 

0.06745 

0.08470 

0.08625 

0.09383 

0.09778 

0.1020 

0.1118 

0.1273 

0.1485 

0.1504 

0.1728 

0.1754 

0.2334 

0.2435 

0.2947 

0.3050 

0.5983 

0.7461 



0.5496 
0.4880 
0.4789 
0.4047 
0.3736 
0.3711 
0.3601 
0.3547 
0.3493 
0.3377 
0,3218 
0.3036 
0.3022 
0.2866 
0.2849 
0.2549 
0.2507 
0.2324 
0.2293 
0.1742 
0.1589 



0.4074 X10- 1 
0.2900X10" 5 
0.2755 X10" S 
0.1 785X10-* 
0.1469x10"' 
0.1447X 10" s 
0.1346X10" 5 
0.1300X10-' 
0.1254X10-' 
0.1159xl0- s 
0.103SX10- 5 
0.9108 X10" 6 
0.901 5 X10" 6 
0.8015X10" 6 
0.7914X10"' 
0.6219 X10- 6 
0.6002 X10- 6 
0.5112X10-" 
0.4967 XIO" 6 
0.282SX10" 6 
0.2354X10-* 



0.4497 
0.5113 
0.5204 
0.5943 
0.6254 
0.6278 
0.6389 
0.6442 
0.6496 
0.6612 
0.6771 
0.6951 
0.6965 
0.7121 
0.7137 
0.7436 
0.7478 
0.7660 
0.7692 
0.8240 
0.8391 



0.6292x10"' 
0.7533 X10- 3 
0.7728 X10" 3 
0.9468 XlO" 3 
0.1029X10" 2 
0.1036X10" 2 
0.1067xl0- 2 
0.1083 XlO" 2 
0.1099X10" 2 
0.1134X10" 2 
0.U83X10" 2 
0.1243X10" 2 
0.1248X10" 2 
0.1303X10" 2 
0.1309xl0- 2 
0.1425X10- 2 
0.1442X10- 2 
0.1 522 XlO" 2 
0.1536X10- 2 
0.1829X10- 2 
0.1929x10-' 



-35.928 
-35.196 
-35.196 
-33.982 
-33.560 
-33.744 
-33.535 
-33.250 
-33.246 
-32.987 
-32.895 
-32.401 
-32.405 
-32.154 
-32.095 
-31.459 
-31.296 
-31.049 
-30.966 
-29.719 
-29.263 



-26.243 
-26.606 
-26.767 
-26.867 
-26.997 
-27.223 
-27.210 
-27.020 
-27.112 
-27.058 
-27.247 
-27.073 
-27.102 
-27.127 
-27.097 
-26.990 
-26.901 
-26.976 
-26.948 
-26.669 
-26.482 



Data of Ramstetter and Hantke [67] 



0.4744 


0.1917 


0.3432X10" 6 


0.S066 


0.1726X10- 2 


-33.6S1 


-30.324 


0.6160 


0.1721 


0.2760X10-" 


0.8260 


0.1843X10" 2 


-33.556 


-30.543 


0.8223 


0.1526 


0.2171X10" 6 


0.8454 


0.1974X10" 2 


-33.388 


-30.718 



Data of Roth and Zeumer [68] 



0.6066 


0.1732 


0.2796X10- 6 


0.8249 


0.1836X10- 2 


-29.522 


-26.490 


0.8566 


0.1500 


0.2099X10" 6 


0,8480 


0.1992X10" 2 


-29.142 


-26.517 


0.9876 


0.1413 


0.1866X10" 6 


0.8566 


0.2058X10- 2 


-28.673 


-26.201 


1.1587 


0.1321 


0.1 635 X10" 6 


0.8658 


0.2133X10" 2 


-28.543 


-26.320 



Data of Johnson and Sunner [37] 



0.01330 


0.6566 


0,7401X10-' 


0.3430 


0.4417x10"' 


-39.212 


-27.623 


0.01416 


0.6466 


0.7003 X10" 5 


0.3529 


0.4579 xlO- 3 


-38.752 


-27.340 


0.01544 


0.6328 


0.6489 XIO" 5 


0.3667 


0.4807X10-' 


-38.376 


-27.209 


0.02272 


0.5708 


0.4626x10"' 


0.4286 


0.5895X10" 3 


-37.074 


-27.013 


0.02932 


0.5301 


0.3706X10-' 


0.4692 


0.6670X10"' 1 


-36.309 


-26.970 


0.03964 


0.4830 


0.2854X10"' 


0.5162 


0.7639X10-' 


-35.384 


-26.882 


0.06039 


0.4204 


0.1986X10"' 


0.5787 


0.9076X10- 5 


-34.342 


-26.950 


0.07829 


0.3842 


0.1590xi0" s 


0.6148 


0.1001 xlO" 2 


-33.731 


-26.98! 


0.09772 


0.3548 


0.I316X10" 5 


0.6441 


O.1083X10- 2 


-33.187 


-26.956 


0.2312 


0.2559 


0.6344x10-* 


0.7427 


0.1421 XlO" 1 


-31.367 


-26.881 
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Figure 2-Experimental enthalpies 
of solution of S0 2 (g) in water at 
298.15 K as a function of 
tx(SO?(aq)). The measurements 
are those of Johnson and Sunner 
[37] ( + ), Stiles and Felsing [36J 
(□), Ramstetter and Hantke [67] 
(0)- and Roth and Zeumer [68] 
(A). The straight line connects 
the selected values of 
A#S=-26.97 and AHi 
= -44.77 k J mol ', 



0.4 0.6 

a (S0 2 (aq)) 

Zambonin and Jordan [38], respectively. Using the 
COD ATA [2] value of 55.8 15 ±0.040 kJ mor 1 for the 
enthalpy of ionization of water, to adjust these re- 
spective values of AH^ leads to -52.8±0.25, 
-50.7+0.32, and -49.5±2.5 kJ mor' for AH§. The 
uncertainties assigned to the measurements of Ram- 
stetter and Hantke [67] and Roth and Zeumer [68] are 
purely statistical. These values of AH? are given in table 
8 together with the value obtained from Zambonin and 
Jordan's [38] measured enthalpy of solution of S0 2 (g) in 



Table 8. Ml at 298.15 K for process (F): S0 2 (g) + H 2 0(/) = 
SO?-(aq) + 2 H+(aq). 



Investigator(s) 



AHVkJ mol" 



Ramstetter and Hantke [67] 
Roth and Zeumer [68] 
Zambonin and Jordan [38] 
Zambonin and Jordan [38] * 
This evaluation 



-52.8±0.23 

-50.7±0.30 

-49.5±2.5 

-47.36±0.9 

-48.42iQ.50 



■ Derived from the authors' measurements of the heat of solution of S0 2 (g) in 
NaHSOj(aq) and the enthalpy data for processes (A) and (B). 



1.C 



dilute NaHS0 3 and their measured values for processes 
(A) and (B). 

It is possible that the result of Ramstetter and Hantke 
[67] are in error by about 4 kJ mol"' since their en- 
thalpies of solution of S0 2 (g) in water (see figs. 1,2) were 
in error by this amount in the same direction as we 
believe these are for process (F), We cannot explain the 
difference of 2.3 kJ mol -1 between the AH° for process 
(F) obtained from Roth and Zeumer [68] and the ten- 
tative selection. As can be seen, the results of Zambonin 
and Jordan [38] bracket this selection. We believe that 
our selected value is best and needs no serious adjust- 
ment. This will become apparent in the next section. A 
modern, accurate measurement of AH" for process (Fl) 
would be of value in confirming our selection. 

12 The Standard State Entropy of SO^-(aq) 
7.2.1 From the Tentative Selections for Process (F) 

Use of the values of AG° tentatively selected for pro- 
cesses (B) and (E) leads to a value of AG£=51.1 1±0.24 
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kJ mol 1 . This value together with AH$ = -48.42 ±0.50 
kJ mol- 1 leads to AS§ = -333. 83± 1.9 J mol" 1 K" 1 . The 
CODATA [2] values for S°(H 2 0,/) (69.950±0.030 J 
mol" 1 K- 1 ) and S (SO 2 ,g) (adjusted to 248.223+0.05 J 
moL 1 K - ' at a standard state pressure of one bar) are 
then used to calculate a value of — 15.66±2.0 J mol -1 
K _1 for the standard state partial molar entropy of 
SOi-(aq). 



7.2.2 From Data on the Na 2 S0 3 System 

There is a direct path to the standard state partial 
molar entropy of S0 3 ~(aq) from the Third Law entropy 
of Na 2 S0 3 (cr) and the value of AG and AH for the 
process: 



Na 2 SOj(cr)= 2 Na + (aq) + S05"(aq). 



(G) 



The stable solid phase in equilibrium with saturated 
Na 2 SOj solution at 298.15 K is the heptahydrate. The 
solubility measurements of Foerster et al. [72] and Kobe 
and Hellwig [73] lead to a solubility of 2.42±0.10 mol 
kg -1 at 298.15 K. The activity of water, 0.908+0.002, 
and the mean ionic activity coefficient, 0.190+0.005, at 
2.42 mol" 1 kg~' are obtained by extrapolation of the 
results of the evaluation of Goldberg [74]. These result 
in AG°=4.016±0.20 kJ mol -1 for the process: 

Na 2 S0 3 -7 H 2 0(cr) = 2 Na + (aq) + SOl-(aq) 

+ 7 H 2 0(/). (Gl) 

Arii [75] reports vapor pressure measurements over 
Na 2 SO r 7 H 2 0(c): 

log p(mm of Hg)=10.656-2797.1/7\ (11) 

These data, after correction for the nonideality of the 
water vapor, and with the CODATA AG° value for the 
vaporization of H 2 0(/), leads to AG°=4.070±0.15 kJ 
mol" 1 for the process: 

Na 2 SO r 7 H 2 0(cr)=Na 2 S0 3 (cr)+7 H 2 O(0- (G2) 

The vapor pressure measurements of Tarassenkow 
[76] over Na 2 SO r 7 H 2 0(cr): 



log p(mm of Hg)= 9.949 -2608/r 



(12) 



were rejected since they lead to an unreasonable AS for 
the decomposition process (G2). Combining AGgi and 
AG& we obtain a value of AG£=-0.054±0.25 kJ 
mol- 1 at 298.15 K. 



Above 308 K the stable solid phase in equilibrium 
with an aqueous solution of Na 2 S0 3 is the anhydrous 
salt. Extrapolation to 298. 1 5 K of the solubility mea- 
surements of Foerster et al. [72] and of Kobe and Hell- 
wig [73] leads to a solubility of 3.26±0.05 mol kg -1 . 
With an exptrapolated [74] mean ionic activity coeffi- 
cient of 0.19±0.03, AG£=0.13±0,6 kJ mor 1 is ob- 
tained, confirming the result obtained from the pathway 
involving the heptahydrate. 

The recent enthalpy of solution measurements of 
Vanderzee [41] yield AHg= - 13.26+0.04 kJ mol -1 and 
are preferred to the earlier measurements of de Forc- 
rand [77] and Kennedy and Lister [78] which yield 
AH£= - 14.9+2.0 kJ mol" 1 and - 18.62±0.4 kJ mor 1 , 
respectively. 

With AHg= - 13.26+0.04 kJ mol- 1 and 
AGS=-0.054±0.25 kJ mol- 1 a value of AS<s = 
-44.29+0.85 kJ mol" 1 is obtained. Use of 
5°{Na 2 S0 3 (cr)} = 145.94±1.2 J moL 1 K" 1 at 298.15 K 
[79] and the CODATA entropy for Na + (aq) of 
58.45+0.15 J mol -1 K _l , results in a partial molar en- 
tropy of -15.25±1.5 J moL 1 K" 1 for SO?"(aq). This 
value is in very good agreement with the value of 
-15.66+2.0 J mor 1 K ' obtained from the S0 2 (g) 
cycle. 



8. The Oxidation of S0 2 (aq) to H 2 S0 4 (aq) 

There are two calorimetric determinations involving 
the oxidation of S0 2 (aq) that were used by the 
CODATA "Key Values" task group [2] as input for 
their evaluation of A f H° of SO^Caq). Reversing their 
procedure, we use the final CODATA value of 
A /J ff a (SOj-(aq))=-909.34±0.40 kJ mor 1 and these 
reactions to obtain AfH"^^^ 2500 H 2 0)). The re- 
actions investigated by Johnson and Sunner [37] and by 
Johnson and Ambrose [80] are, respectively: 

Br 2 (/)+2 H 2 0(/)+S0 2 (in 2500 H 2 0)=H 2 S0 4 
(in 2500 H 2 0) + 2 HBr(in 1250 H 3 0); 
AH=~ 232.09 ±0.42 kJ mol" 1 

and 

Cl 2 (g)+2 H 2 0(/)+S0 2 (in 2500 H 2 0)=H 2 S0 4 
(in 2500 H 2 0) + 2 HCl(in 1250 H 2 0); AH 
= -323.34±0.60kJmor 1 

Using A f H° (HBr(in 1250 H 2 0))= - 121.06+0.15 kJ 
mol- 1 and A f H° (HCl(in 1250 H 2 O))=-166.695±0.10 
kJ moL' from CODATA [2] and Parker [66] and A f H a 
(H 2 S0 4 (in 2500 H 2 O))=-895.43±0.40 kJ mol-' from 
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CODATA and Wu and Young [69] we obtain 
A/if^SOiOn 2500 H 2 0))= -333.80±0.60 kJ mol" 1 and 
-333.82±0.65 kJ mol"', respectively. 

Independently, from the present evaluation we obtain 
A S0] H = - 37. 14±0.20 kJ mol" ' for S0 2 (g)-*S0 2 (in 2500 
H 2 0) which leads to A / £T (SO 2 (in 2500 H 2 0)) = 
— 333.95±0.28 kJ mol -1 . This agreement, well within 
the assigned uncertainties, substantiates the selections 
made here. More importantly, however, the enthalpy 
relationships between aqueous solutions of S0 2 and 
H 2 S0 4 (aq) over a range of concentrations are well de- 
fined. 



9. Final Selected Values for the 
Processes and Properties 

The very near agreement in the partial molar entropy 
of SOj _ (aq) obtained via two independent thermo- 



chemical pathways may be fortuitous and the result of a 
cancellation of errors. Nevertheless, it serves to confirm 
the selections made for the processes used in obtaining 
the partial molar entropy of SOj - (aq). Tables 9 and 10 
contain, respectively, our recommended values for the 
processes and for the thermodynamic properties of the 
species. These values are consistent with the forth- 
coming CODATA "Key Values" [2] adjusted to a stan- 
dard state pressure of one bar. Note that we have adop- 
ted a final "best" value of - 15.40 ±0.8 J mol ~ l K" 1 for 
5°(S03 _ ,aq) and adjusted the value of K B to be 
(6.7 2 ±0.5)xl0- 8 mol kg -1 rather than the value of 
(6.5±0.5)X10 -8 mol kg -1 used in section 4.2. In tables 
9 and 10 the uncertainties have also been adjusted in 
light of the agreement of the data. The thermodynamic 
parameters for process (C) are not as well known as 
those for the other processes. Consequently, the prop- 
erty values for S 2 05~(aq) have much larger uncer- 
tainties than have been assigned to the other species. 



Table 9, Recommended values for the processes involving sulfur(IV)-oxygen species at 298.15 K. 

The value of the equilibrium constant for a given process can be calculated using the equation: R In K~ 

-AG? 98 . l; /298.15 + A/f 2 ° !ls , l3 [l/298.15-l/r] + AC=[(298.15/7')-H-/B(r/298.15)] + (298.15/2)(dAC°/d7')[(?'/298.15)-(298.15/r)-2 
/«( 77298. 15)]. 



Process 



AG 



Aff° 



kJ mol 



AC° 
J mol" 1 K" 



Jmol-'K" 3 



S0 2 (aq) + H 2 O(0 = H + (aq) ± HSOl(aq) (A) 

HSOj(aq) = H + (aq) + SOr(aq) (B) 

2 HSOi(aq) = S 2 Ot-(aq) + H 2 0(/) (C)" 

S0 2 (g) = S0 2 (aq) (D) 

SOj(g) + H 2 0(/) = H + (aq) + HS07(aq) (E) 



10.60+0.10 -17.80+0.40 
40.94±0.20 -3.65±0.10 

8.53±0.80 -4.6±4.0 

-0.51 ±0.10 -26.97+0.30 

10.09+0.14 -44.77±0.40 



-272.+ 10 
-262.±14 
^21.±25 a 
155.++0 
-117.±14 



1.7" 
-2.T 
-1.9 a 
-0.035 3 

+ 1.7" 



J Estimated. 

" For process (C), all values refer to one mole of S,Os ~(aq). 



Table 10. Recommended thermodynamic property values at 298.15 K in SI units and at a standard state pressure of one bar (0.1 MPa). 



Species 


A f H° A / G" 
kJ mol- 1 


5° c; 

J mot"' K"' 


S(cr, rhombic) 
SO a (g) 
SO§(aq) 
S 2 2 5 -(aq) 



-296.81+0.20 8 
-323.78±0.32 
-973.6±4.1 



-300.09±0.21 s 
-3O0.60±0.23 
-S08.61±0.94 


32.054±0.050 a 
248.223 ±0.05" 
159.48±0.75 
154. ±13. 


22.686±0.050 
39.842+0.020 
195. ±10 
-100.0+25" 


HSOl(aq) 
H 2 SO§(aq) e 
SOl-(aq) 
Na 2 S0 3 (cr) 


-627.41 ±0.32 

-609.61 ±0.32 

-631.06+0.40 

-1098.48±0.42 


-527.14±0.25 

-537.74+0.23 

-486.20±0.33 

-1010.08±0.46 

Auxiliary Values" 


134.17±0.65 

229.43+0.75 

-15.40±0.80 

145.94±0.80 


-2. ±10 
270. ±10 
-264. ±10 
120.25 


0»(g) 

H 2 (g) 
H 2 0(/) 

Na(cr) 
Na+(aq) 





-285.830±0.04 


-240.34±0.06 





-237.141+0.04 


-261.95±0.10 


205.152+0.005 

130.680±0.005 

69.950±0.030 

5 1.30 ±0.20 

58.45±0.15 


29.378±0.003 
28.836±0.002 
75.300 
28.23 ±0.20 



' CODATA [2] selections. 

b Estimated. 

c Convention, property values are set equal to the sum of those of SO;Uq) and H;Of/). 
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The property values given in table 10 for HSOl(aq), 
S05~(aq), and S^Oi'Caq) are based upon a more current 
data base than are the property values given in earlier 
evaluations [1,32,81,82]; for this reason they are to be 
preferred. 

Table 11 gives properties (y ± , tj>, L$) and/?(S0 2 ), the 
partial pressure of S0 2 (g) over aqueous solutions of S0 2 
as a function of the stoichiometric molality. The values 
of L$ are obtained as a composite of experimental heat of 
dilution data [40], heat of solution data [36,37], and the 
use of the equilibrium model (sec, 6) to calculate values 
of L$ using the selected property values for processes 
(A) (B), and (C). The properties y±, <p, and L^ all have 
contributions from the various species in the solution. L$ 
is a bulk, stoichiometric property and is the difference in 
enthalpy between the real solution and the hypothetical 
standard state solution consisting of H + (aq) and 
HSOl(aq). The activity and osmotic coefficients were 
calculated using a speciation model of the solution and 
are expressed in table 1 1 for the final treatment of aque- 
ous S0 2 solutions as a uni-univaient electrolyte. The 
excess Gibbs energy is given by G M =v[m s, (H2S03)]RT 



(/« y± — cf> + 1). where v—2. 

Values of y ± can also be calculated from the spectros- 
copic data of Huss and Eckert [24] if one assumes a value 
of K K and a value of the B parameter in eq (5) and applies 
the equation: 



*A=[7(H + )y(HSOy)]m(H + )m(HS07)] 

+ [m(S0 2 )a w ] 



(13) 



to calculate the product [y(H + )y(HSOl)]. Since yf is 
equal to the quantity yfiijmf, the (stoichiometric) mean 
ionic activity coefficient (y+) can be calculated using: 



y ± = [y(H + )y(HSOl)m(H + )m(HSOT)] 
-Kw st (H + )m st (HSOT)] 



(14) 



The average deviation between the values of y ± calcu- 
lated in this way using the spectroscopic data of Huss 
and Eckert [24] and the values of y ± calculated from the 
equilibrium model is 0.0026. 



Table 11. Stoichiometric thermodynamic properties of aqueous sulfur dioxide solutions at 298. 15 K. The values of y -± 
to (H+) (HSOT). 



<fi, and Lt, are relative 



m s, (S0 2 ) 


m 5, (H 2 SOj) 


7= 


* 


L>$ 


/>(S0 2 ) 


mol kg 


-i 






kJ mol"' 
1.09 


bar (0.1 MPa) 


0.001000 


0.001000 


0.909 


0.960 


0.0000483 


0.002000 


0.002000 


0.855 


0.934 


1.91 


0,000171 


0.003000 


0.003000 


0.812 


0.914 


2.58 


0.000347 


0.004000 


0.004000 


0.777 


0.897 


3.14 


0.000566 


0.005000 


0.005000 


0.748 


0.883 


3.61 


0.000818 


0.006000 


0.006001 


0.722 


0.870 


4.05 


0.00110 


0.007000 


0.007001 


0.699 


0.859 


4.43 


0.00140 


0.008000 


0.008001 


0.679 


0.849 


4.78 


0.00173 


0.009000 


0.009001 


0.661 


0.840 


5.08 


0.00207 


0.01000 


0.01000 


0.644 


0.832 


5.36 


0.00243 


0.02000 


0.02000 


0.533 


0.777 


7.36 


0.00664 


0.03000 


0.03001 


0.468 


0.745 


8.54 


0.0115 


0.04000 


0.04002 


0.424 


0.723 


9.36 


0.0168 


0.05000 


0.05004 


0.391 


0.706 


9.99 


0.0224 


0.06000 


0.06006 


0.366 


0.693 


10.46 


0.0282 


0.07000 


0,07008 


0.345 


0.682 


10.85 


0.0341 


0.08000 


0,08011 


0.327 


0.674 


11.17 


0.0402 


0.09000 


0.09014 


0.312 


0.666 


11.44 


0.0464 


0.1000 


0.1001 


0.300 


0.659 


11.66 


0.0527 


0.1500 


0.1504 


0.253 


0.635 


12.56 


0.0850 


0.2000 


0.2007 


0.224 


0.620 


13.14 


0.119 


0.3000 


0.3016 


0.188 


0.600 


13.89 


0.188 


0.4000 


0.4029 


0.165 


0.588 


14.33 


0.258 


0.5000 


0.5045 


0.149 


0.579 


14.64 


0.330 


0.6000 


0.6065 


0.137 


0.572 


14.89 


0.403 


0.7000 


0.7089 


0.127 


0.566 


15.08 


0.477 


0.8000 


0.8116 


0,120 


0.561 


15.24 


0.551 


0.9000 


0,9148 


0.113 


0.556 


15.37 


0.625 


1.0000 


1.0183 


0.107 


0.553 


15.48 


0.700 
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10. The Effects of Perturbations of 

Parameters on Calculated Quantities 

in the Equilibrium Modelling 

0fSO 2 + H 2 O 

The model used to describe the S0 2 + H 2 system 
used several thermodynamic parameters: AG , AH , 
and AC° at 298. 15 K for processes (A), <B), and (C), and 
a B parameter (see eq (4)). It is of interest to examine the 
effects of perturbations in these parameters on the quan- 
tities which can be calculated from the model, namely 
afSO;), y±, tf>, and L*. The examination of the effects of 
these perturbations will be confined to the reference 
temperature of 298.15 K. 

The results of these calculations are summarized in 
table 12, Each parameter has been perturbed by the limit 
of error assigned to that parameter in table 9; the B 
parameter has been perturbed by 0.5. As a result of these 
calculations, the values of K A and B are found to be most 
significant in calculating value of a(SO?), y±, and d>. In 



the calculation of L$, the most important quantity is 
AHt followed by K K and B, 

The results of these calculations are useful in esti- 
mating uncertainties in the thermodynamic properties 
that were calculated from the equilibrium model. Thus 
the calculated values of K D are uncertain by at least 1 to 
3% because of uncertainties in the calculation of a(SO§). 
The uncertainties in y±, §, and £,$ are given in the last 
column of table 12; uncertainties in y± range from 0.2 to 
1.4%, in <f> from 0.1 to 0.3%, and in L$ from 2 to 3%. 
This sensitivity analysis does not consider possible cou- 
pling effects between the parameters varied. 



11. Extension of the Temperature Range 

Over Which Properties and 

Equilibrium can be Calculated 

In the model used, it was assumed that the heat capac- 
ities of the species were independent of temperature. 
This was necessary since we have no direct knowledge 
of the temperature dependency of the heat capacities of 



Table 12. Percentage effects on calculated quantities (ct(SO§), -y*, 4>, andi*) due to perturbations in the parameters of the model 
used to describe the thermodynamics of aqueous SOi solutions. The parameters which were perturbed were B in eq (4) (1.5 instead of 2.0), 
K A (0.0143 instead of 0.0139 mol kg" 1 ), -Kb (7.2 X 10 -8 mol kg" 1 instead of 6.7 X 10" s mol kg -1 ), K c (0.042 instead of 0.032), AHJL (- 18,20 
instead of - 17.80 kj mol" 1 ). A#S (-3.75 instead of -3.65 kJ mo)-'), and AH£ (-8.60 instead of -4.60 kJ mol" 1 ). The temperature is 
298.15 K in all cases, A discussion of the effects of the B parameter on the value of A&Z is given in section 4.1. 



m 








Parameters Modified 








total in 


mol kg -1 


B 


K K 


Kb 


K c 


AH% 


Am 


AHZ 


quadrature 










Effects on a(SO§) 










0.001 


-0.08 


-2.5 


-3.3X10- S 


-0.0018 











2.5 


0.01 


-0,34 


-1.5 


-8.3X10-" 


-0.0093 


— 


— 


— 


1.5 


0.10 


-0.50 


-0.63 


-1.3XIO" 6 


-0.022 


— 


— 


— 


0.80 


1.00 


-0.47 


-0.23 


-2.0X10-' 


-0.035 
Effects on 7* 








0.52 


0.001 


-0.044 


+0.16 


-1.7X10- 5 


-0.00091 











0.16 


0.01 


-0.17 


+0.64 


-3.8X10" 6 


-0.0047 


— 


— 


— 


0.66 


0,10 


-0.25 


+ 1.1 


-5.8X10- 1 


-0.011 


— 


— 


— 


1.1 


1.00 


-0.23 


+ 1.3 


-8.5X10- S 


-0.017 
Effects on d> 








1.3 


0.001 


-0.021 


+0.075 


+ 0.028 


-0.00045 











0.08 


0.01 


-0.07 


+0.26 


+ 3.5X10- 5 


-0.00045 


— 


— 


— 


0.27 


0.10 


-0.05 


+0.28 


+ 6.0X10" 6 


-0.0040 


— 


— 


— 


0.28 


1.00 


-0.021 


-0.15 


+ 1.0XKT 6 


-0.0050 
Effects on Lj, 








0.15 


0.001 


-0.028 


-2.4 


-I.8X10" 3 


-0.0056 


+2.1 


-6.6X10- 4 


-0.011 


3.2 


0.01 


-0,27 


-1.5 


-4.8X10- 5 


-0.014 


+2.3 


-1.6X10' 5 


-0.014 


2.8 


0.10 


-0.44 


-0.61 


-3.5X10" 6 


-0,027 


+ 2.3 


-l.ixio-* 


-0.018 


2.4 


1,00 


-0.42 


-0.22 


-5.0X10" 7 


-0.041 


+ 2,3 


-1.3X10- 7 


-0.025 


2.4 
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these species. However, using the "correspondence 
principle" of Criss and Cobble [83] and the parameters 
given by them in their table 1, the quantitiy (dC°/dT) 
can be estimated for the aqueous species HSOj, SO 2 ~, 
and S 2 Oi~. Doing this, we obtain (dC;/dr)= + 1.7, 
-1.0, and +1.5 J mo!" 1 K 2 for HSO7, SO 2 ,", and 
S 2 2_ , respectively. From data on the heat capacity of 
water and of S0 2 (g) [7,82], {dC°/dT)= +0.0066 J 
mol" 1 K- 2 for HjO(/) and (dC°/dT) = +0.035 J mol -1 
K- 2 for S0 2 (g); we estimate (dC;/dF) = for SO?(aq). 
These estimates are combined to yield values of +1.7, 
-2.7, -1.9, -0.035, and +1.7 J mol 1 K" 2 for 
{dAC°/d7) for processes (A), (B), (C), (D), and (E), 
respectively. These values are also given in table 9. 

If eq (6) is used without including a (dACp/dT) term, 
values of AT A =0.00153 mol kg -1 and i<: D =0.213 3 mol 
kg"' bar -1 at 373.15 K are calculated. Inclusion of the 
(dACp/d T) terms and the estimates of them given above 
leads to K A =0.0017 1 mol kg" 1 and K D = 0.2 12„ mol kg" 1 
bar -1 at 373.15 K. Thus, over the temperature range 273 



to 373 K the thermodynamics of the S0 2 + H 2 system 
are reasonably well described by the parameters given 
in tables 9 and 10. A plot of the measured Henry's Law 
constants as a function of temperature and the curve 
calculated using the final selected values given in table 
9 are shown in figure 3. A particularly useful series of 
experiments would be the measurement of the heat ca- 
pacities of aqueous sulfur dioxide solutions from 273 K 
to temperatures greater than 373 K. 



We thank Drs. Celina Dobrogowska, Loren Hepler, 
Peter Tremaine, and Cecil Vanderzee for sharing the 
results of their research with us prior to publication and 
for their helpful discussions. 




Figure 3-The measured Henry's 
Law constants as a function of 
temperature and the curve (solid 
line) calculated using the final 
selected values in table 9. The 
data sets and their correspond- 
ing symbols are: (.0) Hudson 
[60]; (4M Maass and Maass [61]; 
(O) Campbell and Maass [29]; 
I©') Morgan and Maass [31]; (Q) 
Conrad and Beuschlein [58]; (*) 
Johnstone and Leppla [13]; (D, 
O, A, and +) Beuschlein and 
Simenson [56J at, respectively, 
0.51, 1.09, 4.36, and 7.45 g 
S0 2 / 1 00 g H 2 0; (ZK) Rabe and 
Harris [62]; (SO Vosolobe et al. 
[64]; (ffl) Tokunaga [63]; (y) 
Douabul and Reilly [59]; and 
(X)Byerley [57]. 
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Glossary 



Roman 



a 

f 
m 
n 

P 

z 

A m 

B 

C P 

G 

H 

I 

K 

L, 

R 
S 

T 



activity 

fugacity 

molality/mol kg -1 

amount or number of moles of substance 

pressure 

charge 

Debye-Huckel constant; A m = 1.17642 kg m 

mol _!/2 at 298.15 K 

parameter in Debye-Hiickel equation 

heat capacity at constant pressure 

Gibbs energy 

enthalpy 

ionic strength 

equilibrium constant 

relative apparent molar enthalpy, equal to 

H"/m 

gas constant; R =8.31448 J mor 1 K~ l 

entropy 

temperature 



Greek 

a fraction 

y activity coefficient 

<p osmotic coefficient 

A change 

Superscripts 

ex excess 

st stoichiometric 

o value of a quantity which applies 

the thermodynamic standard state 
° designates a neutral species, i.e., SO? 

species quantity 



at 



Subscripts 

/ formation property 

4> apparent molar quantity 

± mean ionic property 

w designates water 
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